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We study the ground state phases of Bose-Fermi mixtures in one- dimensional optical lattices 
with quantum Monte Carlo simulations using the Canonical Worm algorithm. Depending on the 
filling of bosons and fermions, and the on-site intra- and inter-species interaction, different kinds of 
incompressible and superfluid phases appear. On the compressible side, correlations between bosons 
and fermions can lead to a distinctive behavior of the bosonic superfluid density and the fermionic 
stiffness, as well as of the equal-time Green functions, which allow one to identify regions where the 
two species exhibit anticorrelated flow. We present here complete phase diagrams for these systems 
at different fillings and as a function of the interaction parameters. 
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I. INTRODUCTION 

The experimental realization of strongly correlated sys- 
tems with ultracold gases loaded in optical lattices [l| has 
generated tremendous excitement during recent years. 
Initially thought of as a way to simulate condensed mat- 
ter model Hamiltonians, like the Bose-Hubbard Hamil- 
tonian 0, loading atoms on optical lattices has enabled 
the creation of quantum systems that are unexpected 
in the condensed matter context. Among these systems 
the realization of Bose-Fermi mixtures in optical lattices 
0j 0i IE] J where the inter- and intra-species interactions 
can be tuned to be attractive or repulsive is a re- 
markable example of the scope of realizable models. 

Theoretical studies of Bose-Fermi mixtures in one- 
dimensional lattices have been done for homo gen eous 
0, S S M M M M, [13 and trapped [ll 111 [3 
systems. Several approaches have been used: Gutzwiller 
mean- field theory [IJ] , strong coupling expansions 
bosonization [7 
ical [nl, [H, \n 



and exact analytical [10[ and numer- 
161 . studies. Recently, a mixture of 
bosonic atoms and molecules on a lattice was studied 
numerically . The landscape of phases encountered is 
expansive, and includes Mott insulators, spin and charge 
density waves, a variety of superfluids, phase separation, 
and Wigner crystals. However, the phase diagram in the 
chemical potential-interaction strength plane has not yet 
been reported. 

It is our goal in this paper to present a study of re- 
pulsive Bose-Fermi mixtures in onc-dimcnsional lattices 
that generalizes previous studies, which focused on spe- 
cific special densities, to more general filling. After map- 
ping the phase diagram we will explore different sec- 
tions in greater detail. Since the particular case in which 
the lattice is half filled with bosons and half filled with 




fermions has been carefully studied in Ref. we will 
instead concentrate here on two cases: (i) when the num- 
ber of bosons is commensurate with the lattice size but 
the number of fermions is not, and (ii) when the sum of 
both species is commensurate with the lattice size but 
the number of bosons and fermions are different. Some 
of the phases present in these cases have been identified 
by Scngupta and Pryadko in their grand canonical study 
in Ref. |12l | and by Hebert et al. in the canonical study 
recently presented in Ref. ■ 

The Hamiltonian of Bose-Fermi mixtures in one dimen- 
sion can be written as 
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where {b^ ) are the boson creation (destruction) opera- 
tors on site I of the one-dimensional lattice with L sites. 
Similarly, //(/; ) are the creation (destruction) operators 
on site I for spinless fermions on the same lattice. For 
these creation and destruction operators nf'^ are the 
associated number operators. The bosonic and fermionic 
hopping parameters are denoted by ts and tp respec- 
tively, and the on-site boson-boson and boson-fermion 
interactions by Ub b and Ub f ■ In this paper we will con- 
sider the case Ib = tp = 1 (i-e. when the boson and 
fermion hopping integrals are equal) and choose = 1 
to set the scale of energy. 

It is useful to begin a discussion of the phase diagram 
with an analysis of the zero hopping limit (Ib = t p = 0) 
similar to the one done by Fisher et al. in Ref. [101 for 
the purely bosonic case. Consider a particular fermion 
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occupation of one fourth of the lattice sites, Np = L/4 
fixed. Bosons can be added up to A''^ = 3L/4 without 
sitting on a site which is aheady occupied by either a 
boson or a fermion. Therefore the associated chemical 
potential /i is small. What happens when Nb exceeds 
3i/4 depends on the relative strength of Ubb and Ubf- 
li Ubf is less than 2Ubb then the extra bosons sit 
atop of the fermions and /i jumps by Ubf- The chem- 
ical potential stays at this elevated value of Ubf until 
all the sites with fermions also have a boson. At that 
point additional bosons start going onto sites with a bo- 
son already, and /i jumps to 2Ubb- Thus in general there 
are incompressible phases where the boson chemical po- 
tential jumps both at commensurate pB = 1, 2, 3, ... (as 
for the pure boson-Hubbard model) and also at = 
1 — pF, '2 — pF, 3 — pf, For [/si? greater than 2L/bb 
and less than 6Ubb the incompressible phases still start 
&t pb ~ ^ — Pf but the following potential jumps are 
shifted up by 1 — pp- Turning on the hoppings tB,tp in- 
troduces quantum fluctuations which will ultimately de- 
stroy these Mott plateaus and introduce new, intricate 
phases. 

II. CANONICAL WORM ALGORITHM 

We perform Quantum Monte Carlo simulations 
(QMC) using a recently proposed Canonical Worm al- 
gorithm [lO, |2l|. This approach makes use of global 
moves to update the configurations, samples the wind- 
ing number, and gives access to the measurement of n- 
body Green functions. It also has the useful property of 
working in the canonical ensemble. This is particularly 
important for the present application since working with 
two species of particles leads to two different chemical 
potentials in the grand canonical ensemble. These prove 
difficult to adjust such that the precise, desired fillings 
are achieved. In our canonical simulations the Bose and 
Fermi occupations are exactly specified and the chemical 
potentials /ib and pp are instead computed via ap- 
propriate numerical derivatives of the resultant ground 
state energy (e.g. pB = Eq{Nb + 1) - Eo{Nb))- 

The Canonical Worm algorithm is a variation of the 
Prokof'ev et al. grand-canonical worm algorithm [2^. 
Within the Canonical Worm approach one starts by writ- 
ing the Hamiltonian as H = V — T, where T is comprised 
of the non-diagonal terms and is by necessity positive def- 
inite. The partition function Z = Tre^''^ takes the form 

Z = Tre-^^T,e/o'*(^)'^^ (2) 
= Tre-''^^ [f{Tn)---f{n)dn---dTn (3) 

„ J0<Ti<-<T„<l3 

where T{t) = e^^Te~'^^ . In order to sample expres- 
sion Q an extended partition function is considered by 
breaking up the propagator at imaginary time r and in- 
troducing a "worm operator" W = ^ijki'^ijkib\bjflfi 



that leads to Z{t) = Tre^^'^^^'^M^e"^^. Complete sets 
of states arc introduced between consecutive T operators 
to allow a mapping of the ID quantum problem onto 
a 2D classical problem where a standard Monte Carlo 
technique can be applied. Measurements can be per- 
formed when configurations resulting in diagonal matrix 
elements of W occur. This way unphysical movements 
are exploited to help explore the Hilbert space, but are 
ignored when sampling for measurements. 

As with pure bosonic systems, the evolution of the 
boson and fermion densities pb,Pf with the associated 
chemical potential pbiI^f identifies Mott insulating be- 
havior [l9|. A jump in p signals a Mott phase where 
the compressibility kb = dpB/dpB or Kp = dpp/dpp 
vanishes. 
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FIG. 1: (Color online) ps as a function of chemical potential 
Pb- The fermion density is pp = 1/4 and the interaction 
strengths are fixed at Ubb = 10 and Ubf = 16. There are 
Mott plateaus at pB = 1 — pp = 3/4 and ps = 1 as predicted 
by the tB = tp = analysis. The positions of the Mott lobes 
coincide for different lattice sizes L = 36, 44 and temperatures 
(3 = 32, 48 to within our error bars, which are smaller than 
the symbol size. The dependence of ps on ps in the absence 
of fermions is given for comparison. 

Quantities of interest that we measure include the 
bosonic superfluid density and the fermionic stiffness, 

Pl ^ (VF|)L/2/3 

p^p = (VF|)L/2/3. (4) 

Here {W'^) are the associated winding numbers. Corre- 
lations between the bosonic and fermionic winding num- 
bers pd| are determined by the combinations, 

Pc = {{WB + Wpf)L/2f3 

Pi = {{WB~Wpf)L/2(3 . (5) 

In addition to the usual bosonic and fermionic Green 
function, 

Gf, = (bib,) 

G^^iflf,), (6) 
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we also measure the composite anti-correlated two-body 
Green function 



ij 



(7) 



In Gfj, the fcrmion and boson propagate in opposite di- 
rections (one from j to i and one from i to j). 

The Fourier transforms of Gfj and give the densi- 
ties nB{k) and np{k) in momentum space; na{k) is the 
Fourier transform of the composite two-body Green func- 
tion Ciy 

We performed extensive checks of the code against 
other quantum Monte Carlo simulations in the pure bo- 
son and pure fcrmion cases, and against exact diagonal- 
ization and Lanczos calculations for mixed systems on 
small lattices. 



III. PHASE DIAGRAM IN THE iib-Ubf PLANE 

We begin our determination of the phase diagram by 
calculating the dependence of the density ps on chemi- 
cal potential /i^, mapping out the extent that the Mott 
plateaus described in the introduction survive the intro- 
duction of quantum fluctuations tg^tp. We examine a 
system with a fixed Ubb = 10 and pp ~ 1/4 and focus 
on the regions through ps < 3/2 and Ubp < 5Ubb/2 in 
the phase diagram. The ts = tp = analysis suggests 
for Ubf < '2Ubb there will be plateaus with compress- 
ibility K = at + = 1 (i-e. pb = 3/4) caused by 
Ubf and at = 1 caused by Ubb- Fig. [T] exhibits these 
plateaus for Ubp ~ 16 and ts ^ tp = 1. The complete 
phase diagram in the pb-Ubf plane at fixed Ubb = 10 
is obtained by replicating Fig. [T] for different UbF: and 
is given in Fig. For weak Ubf the phase diagram 
is dominated by the = 1 plateau where the chemical 
potential jumps by 2Ubb — Ubf ~ 2[/bb = 20. As Ubp 
increases, this plateau shrinks and finally terminates at 
Ubf ~ "^Ubb = 20. At the same time, the plateau at 
Pb = I ~ Pf grows to Ubf = 20. The explanation of the 
labeling of the different phases (I- VI) will be given after 
we discuss the supcrfluid response of the system. Fig. [2)3 
shows the phase diagram in the pB and Ubb plane. 



IV. SUPERFLUID RESPONSE AT ps + pf = 1 

After determining the positions of the Mott plateaus, 
we examine the stiffness and Green functions. We take a 
"horizontal" cut through Fig. [5^ by fixing pB + pp = 1 
{pb = 3/4) and increasing Ubp- In Fig. [3] we see that 
for Ubf ^ 2Ubb = 20 the interaction strength Ubf is 
small enough that fermions and bosons can briefly in- 
habit the same site. Now, when a boson visits the site of 
a neighboring fcrmion (or vice versa) it is equally likely 
that the fcrmion will exchange as for the boson to return 
to its original site. Through these exchanges the bosons 
and fermions can achieve anti-correlated winding around 
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FIG. 2: (Color online) Top panel: Phase diagram in the ps- 
Ubf plane for Ubb = 10 and pp = 1/4 obtained by a se- 
quence of plots such as that in Fig. [1] The vertical line at 
Ubf = 16 corresponds to the coupling value in Fig.[T] Within 
the regions labeled pB = 0.75 and ps = 1, the boson density is 
frozen even though the chemical potential varies. The ps = 1 
lobe is pinched off at Ubf ~ 2Ubb- The labeling of the 
phases is: I. Superfluid (ps 7^ p^ 7^ 0); II. Anti- Correlated 
phase (p^g 7^ 0, p^ 7^ 0, Ps = Pf); HI- Anti- Correlated phase 
/ Relay Superfluid (p'^ 7^ 0, p|? 7^ 0) (see text for detailed 
explanation); IV. Mott Insulator / Luttinger liquid {p% — 0, 
Pf 7^ 0)i V. Insulator, ps = pi? = 0; VI. Phase separation. 
Bottom panel: Phase diagram in the pb-Ub b plane for Ubf = 
10 and pF = 1/4. The upper boundary of phase // is de- 
fined by the increase of energy when a boson is added at 
filling Pb + Pf = 1; this is approximately ps = Ubf (for 
Ubf < 2Ubb), a line of slope 1 in the pb-Ubf plane and 
slope in the pb-Ubb plane. A similar strong coupling anal- 
ysis applies for the other boundaries. (See also Fig. [13] and 



the lattice. Thus, the bosonic superfluid density and the 
fermionic stiffness are both non-zero and identical [25| . 
However, as Ubf increases past Ubp ^ 2Ubb = 20 the 
cost of double occupancy becomes prohibitive. With its 
benefits outweighed by energy penalties exacted by Ub f , 
all anti-correlated "superfiuidity" ceases. PoUet et ai 
[Tl] | have argued that this region exhibits phase separa- 
tion. Indeed we do detect a signal of phase separation 
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FIG. 3: (Color online) "Horizontal" sweep, that is fixed den- 
sities pb ~ 3/4 and pp = 1/4 and varying Ubf, through the 
phase diagram of Fig. [2] As in Figs. [T] and O Ubb — 10. 
Top panel: Both the bosonic and fermionic species exhibit 
a finite stiffness at weak coupling (region II of phase dia- 
gram), which decays as Ubf increases until insulating be- 
havior occurs (region V of phase diagram). Bottom panel: 
Near Ubf ~ both pi and p^ are very similar. 
However, quickly after turning on Ubf, the correlated and 
anti-correlated stiffness show that the bosons and fermions 
propagate in opposite directions 7^ 0, while p'l — 0. 



through density structure factor. But the signal is weak, 
about 20 times weaker than what we get at phase VI 
(next section), and compressibility is close to zero, so we 
label this region as an insulator. 

From the results depicted in Fig. [3] one should notice 
that while for quantities like the energy and Mott gap 
/? = 32 is sufficiently low for L = 32 to capture the 
ground state behavior, for stiffnesses one requires much 
lower temperatures. 



V. SUPERFLUID RESPONSE AT ps = 1 

Although it shares the property that kb — with the 
Pb + Pf = ^ lobe, a 'horizontal' cut (Fig. 4) through the 
Pb = I Mott lobe exhibits rather different superfluid re- 




FIG. 4: (Color online) Same "horizontal cut" as Fig.|3]except 
at commensurate density for the boson species alone ps = 1 
and Pf = 1/4. Top panel: The bosons are insulating at weak 
Ubf within this Mott lobe of commensurate bosonic filling 
(region IV of phase diagram) . However the fermions are free 
to flow on the uniform boson background and have nonzero 
stiffness. Upon emerging from the lobe, at Ubf ~ 2Ubb, 
p% becomes nonzero in a window where the two repulsions 
work against each other. After the peak, the system becomes 
phase separated and p^B ^-nd pj- go to zero. Bottom panel: The 
correlated and anti-correlated stiffnesses are essentially equal 
throughout the weak coupling because the flow is dominated 
by fermions. However, in the window the anti-correlated stiff- 
ness increases beyond the correlated stiffness in a weak simu- 
lacrum of phase II (as discussed in the text). 



spouse. This trajectory initially lies within the Mott lobe 
and then emerges into a region of non-zero compressibili- 
ties. As expected, the plateau in pB (Mott gap) indicates 
the bosons are locked into place by the strong Ub b , and 
as a consequence = (Fig. |4]). Throughout this bo- 
son Mott lobe the fermions are, however, free to slide 
over the bosons and so p|, is non-zero. In this region, as 
expected, the fermion compressibility Kp is nonzero. 

The Bose- Fermi repulsion Ubf competes with Ubb 
and, in a window around Ubf ~ '^Ubb, it is energeti- 
cally equivalent for a boson to share a site with another 
boson as with a fermion. The Mott lobe is terminated 
and a superfluid window opens for both species. Finally, 
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for Ubf > '^Ubb, it is energetically unfavorable for a bo- 
son to share a site with a fermion. Wc enter a region of 
phase separation where superflow for both species stops, 
but the compressibilities kb and Kp are nonzero. We 
also confirm phase separation through a density struc- 
ture factor. See also [isLflTj. 



VI. SUPERFLUID RESPONSE AT GENERAL 
FILLING 

Further insight into the physics of this phase diagram 
can be obtained by measuring the superfluid response 
along the same 'vertical' cuts through the phase diagram 
as done in Figs. [T] and [21 in which is varied at fixed 
Ubf- In Figs. [5] and [SI we show the result. Distinctive 
densities in the latter figures are pB — 3/4 (so that pB + 
Pf = 1) and pb = 1. We discuss first (Fig. [S]) the case 
of Ubf = 16, where increasing pB cuts through both 
Mott lobes. The bosonic superfluid density vanishes at 
Pb — 1; dips at Pb + Pf ~ smd is nonzero above, 
below, and between the lobes. The fermion superfluid 
density is never driven to zero in this cut, and only dips 
at the special value pb = 3/4 where the commensurate 
total density works against superfluidity. In the case of 
Ubf = 24, Fig.[6l as pB increases we cut through only the 
Pb = 3/4 lobe. Here the superfluid density is pushed to 
zero for the entire region between = 3/4 and = 1, 
and is non-zero without. 

We now fill in the labeling of the phase diagram of 
Fig. [21 A gapless superfluid phase (I) with (p^ ^ and 
Pp ^ 0) exists at low fllhng of the lattice pb + Pf < 1- 
When the combined fllling of the two species becomes 
commensurate, an anti-corrclated (II) phase appears in 
which p^g and Pp ^ 0, but pg = p^p. This phase 
is characterized by superflow of the two species in oppo- 
site directions and is gapped to the addition of bosons or 
fermions. The usual bosonic Mott insulator, phase IV, 
occurs at commensurate boson densities. However, it can 
be melted by increasing Ubf since the jump in bosonic 
chemical potential (Mott gap) is reduced to 2Ubb — Ubf- 
There is no jump in iip. Eventually quantum fluctu- 
ations break this gap and superflow is allowed. When 
Ubf exceeds 2Ubb, all superflow stops and we enter the 
insulating region V of the phase diagram. 

We speculate that the nature of the superfluidity in 
the narrow phase III, which exists between the two Mott 
lobes is an unusual "relay" process. It is similar to the 
usual superfluid which exists between Mott lobes in the 
single species model, in that p% ^ 0. However, the 
temperature scale at which superfluid correlations build 
up is dramatically reduced. This occurs because the 
bosons can exhibit superflow only by traveling along with 
a fermion partner, and being handed off from fermion 
to fermion in order to wind around the entire lattice. 
The point is that because 2Ub b exceeds Ub f the bosons 
doped into the lattice above pB = I — pp = ■i/4: are forced 
to sit on a fermion. They cannot hop off, but the fermion 
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FIG. 5: (Color online) "Vertical" sweep across ps at pp = 
1/4, Ubb = 10 and Ubf = 16. Top panel: The boson super- 
fluid density changes much less as (3 is increased from /3 — 32 
to /3 = 108 in the superfluid phase I at p_B < 3/4 and pB > 1 
than for phase III 3/4 < pB < 1. This is a hallmark of the 
"relay" superfluid discussed in the text. Bottom panel: The 
correlated winding decreases to zero at ps ~ 3/4 while the 
anti-correlated winding remains flnite. As pB is increased be- 
yond 3/4 the correlated winding increases and overtakes the 
anti-correlated winding. This is another sign of the "relay" 
superfluid. 



can move since it has already paid Ubf share a site 
with a boson. Now, the fermions cannot pass each other 
once a fermion riding atop bosons runs into a fermion 
alone on a site. The fermion without a boson cannot 
move out of the other fermion's way either. However, the 
boson sharing a site with the mobile fermion can then 
hop to the immobile fermion at no energy cost. Thus, 
the boson is passed from one fermion to the other, grant- 
ing it mobility. Signatures of this phase are the lower 
value of the temperature at which the superfluid density 
builds up, that > pg, and more correlated winding 
than anti-correlated. However, there is nothing prevent- 
ing lone fermions from acting as in the anti-correlated 
superfluid phase. Unfortunately this means that poten- 
tial signals are masked. While we do see some of these 
signatures (Fig. [S]) in the specified region, the numbers 
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FIG. 6: (Color online) "Vertical" sweep across ps at pp = 
1/4, Ubb = 10 and Ubf ~ 24. Unlike the case for the weaker 
coupling Ubf ~ 16 in Fig. O both superfluid densities vanish 
in the insulating region of the Mott lobe at commensurate 
total filling. 



are not completely conclusive and will require further in- 
vestigation. 



VII. MOMENTUM DISTRIBUTION 
FUNCTIONS 

To further explore the nature of the phases we turn to 
the momentum distributions for the bosons, fermions, 
and anti-correlated pairing - Fig. [7] - [T^ Each plot 
is made at P — 108 and correspond to the parameter 
choices: I. Ubf ~ 16 and Ng = 20; II. Ubf = 16 and 
Nb = 27; III. Ubf = 16 and Nb = 32; IV. Ubf = 16 
and Nb = 36; V.Ubf = 30 and Nb = 27; VI. Ubf = 28 
and Nb = 36. In the superfluid phase (I.), there is 
a peak in the boson distribution and a plateau in the 
fermion distribution, implying quasi-condensation in the 
bosonic sector and Luttinger liquid-like behavior in the 
fermionic one. In the anti-correlated phase (II.) there is 
neither of the former behaviors, but the Fourier trans- 
form of the anti-correlated pairing Green function has 
a clear Fermi momentum showing Luttinger-like physics 




FIG. 7: (Color online) I. Ubf = 16, Ubb = W, P = 108 and 
Nb ~ 20; Momentum distributions for bosons and fermions, 
and Fourier transform of the anti-correlated pairing (na(fc)) 
Green function. The sharp peak in bosonic momentum dis- 
tribution indicates the presence of a quasi-condensate, while 
fermions have a plateau indicating Luttinger liquid like be- 
havior with a clear Fermi momentum, a property that is 
also shared by the composite fermions described by the anti- 
correlated pairing. 
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FIG. 8: (Color online) II. Ubf = 16, Ubb = W, (3 = 108 
and Nb ~ 27. Bosons do not have a peak at fc = and the 
Fermi momentum is washed out, both reflecting the onset of 
short range one-particle correlations. On the other hand, the 
plateau in the Fourier transform of the anticorrelated pair- 
ing shows that the composite fermions formed by pairing a 
fermion and a boson have a well defined Fermi momentum 

El. 



of the composite fermions (formed by pairing a fermions 
and a boson) [lll,!!^. The "relay" superfluid phase (III.) 
displays momentum distributions that are similar to the 
ones of superfluid phase (I). Next, in the Mott insula- 
tor / Luttinger liquid phase (IV.) one can see a clear 
Fermi momentum in the bare fermion nF^k) and a very 
smooth behavior of nsik) and na(fc), which show that 
their real space Green function counterparts are decay- 
ing exponentially. In the insulating phase (V.) all the 
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FIG. 9: (Color online) III. Ubf = 16, Ubb = W, f3 ^ 108 
and Nb = 32. Qualitatively, this picture is similar to Fig. [7] 
- we have a peak in the bosonic momentum distribution and 
a plateau in the fermionic and anticorrelated pairing momen- 
tum distributions, all indicating power-law decaying correla- 
tions of their corresponding real space Green functions. 
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FIG. 11: (Color online) V. Ubf = 30, Ubb = 10, /3 = 108 
and Nb = 27. The momentum distribution functions in this 
case exhibit the behavior expected from an insulator, i.e., no 
sharp peak in UBik), no plateau in npik), and no Fermi edge 
in na{k). 
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FIG. 10: (Color online) IV. Ubf ^ 16, Ubb = W, P ^ 108 
and Nb = 36. There is no sharp peak in nB{k) and a plateau 
in npik) is present. This phase is a Mott Insulator for bosons 
and Luttinger liquid behavior for the fermions. In this case 
the composite fermions do not exhibit a Fermi momentum. 



correlations decay exponentially and their corresponding 
momentum distribution functions arc smooth functions 
of k. In the case of phase separation (VI.) the bosonic 
momentum distribution is similar to the superfluid, while 
fermionic distribution is insulating. 



VIII. CONNECTION TO PREVIOUS 
THEORETICAL WORK 

As reviewed in the introduction, there is an extensive 
theoretical literature on Bose- Fermi mixtures. We now 
make more detailed contact with previous work, first by 
comparing our results to the strong coupling phase di- 
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FIG. 12: (Color online) VI. Ubf = 28, Ubb = 10, /3 = 108 
and Nb ~ 36. Phase separation. The boson momentum 
distribution function has a peak indicating that there may be 
a kind of superflow in their separate area. Fermions, on the 
other hand, behave as an isolator. The na{k) curve indicates 
that the coupling between bosons and fermions is weak, as we 
would expect in phase separation. 



agram of Lcwcnstcin et al. (LSBF) Q. Fig. 1131 com- 
bines our results and those of LSBF. Besides the quanti- 
tative agreement, we note the following correspondences: 
LSBF's region < /2 < 1 is analogous to our < /i < 20, 
and < a < 1 to our Q < Ubf ^ 20. Furthermore, 
our phase // (Anti-Correlated phase) corresponds to 
LSBF's phase II fl (Fermi liquid of composite fermions 
formed by one bare fermion and bosonic hole); our phase 
IV (Mott insulator / Luttinger liquid) to LSBF's phase 
Ifl (Fermi liquid); and finally our phase /// (Anti- 
Correlated phase / Relay superfluid) to LSBF's phase 
Idw (Density wave phase). These three phases have 
similar qualities and occur approximately at the same 
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FIG. 13: (Color online) A comparison of our phase diagram 
(Fig. ^ with the strong coupling boundaries. Symbols and 
dashed lines are the results of the present QMC work, while 
the solid lines are for ~ tp ~ 0- Unsubscripted Roman 
symbols denote our phases while subscripted Roman symbols 
are the labeling of Lewenstein et al. §]. 

locations at our and LSBF's phase diagrams. 

Both calculations suggest the existence of compos- 
ite particles. Our phase V (Insulator) corresponds to 
LSBF's phase II fd, a region of fcrmionic domains of 
composite fermions formed by one bare fermion and 
bosonic hole. There is one case when the phases do 
not seem to correspond well, namely LSBF's phase II sf 
which is a superfluid of composite fermions formed by 
one bare fermion and bosonic hole. Our results (Fig. 3) 
instead suggest that in this region of Ubf > '2Ubb, the 
superfluid densities vanish, or arc very small. 

IX. EXPERIMENTAL ISSUES 

Albus et al. have given the correspondence 

between Hubbard model parameters Ub B,UBF,tB,tF 
and experimentally controlled parameters. Ubb and 
Ubf are determined by the optical lattice depth, laser 
wavelength, and harmonic oscillator lengths, as well as 
by the scattering lengths gbb and qbf which can be 



tuned by traversing a Feshbach resonance. Similarly, the 
hoppings Ib and follow from the lattice depth and 
atomic masses. It is possible to choose experimentally 
reasonable values of these parameters to correspond to 
the energy scales chosen in our paper. For example, 
following Albus et al, for a ^"^ Rh,'^^ K mixture and laser 
wavelength 600rtm, ubb ~ lOOao, aBF = 123.74ao, 
and Vb = 0.7614 in units of boson recoil energy, with 
Ig = 17.04717TI, we get in units of i^: tB = 1, = 2, 
and Ubb = Ubf = 10. In this paper we have used 
tB ~ tp = 1 , which would be accessible in mixtures with 
mF ~ rriB such as '^'^K,'^^ K. 

X. CONCLUSIONS 



In conclusion, we have mapped out the boson density 
- interaction strength phase diagram of Bose-Fermi mix- 
tures. The Mott lobe at commensurate total density has 
nontrivial superfluid properties, where the two compo- 
nents of superflow can be nonzero and anti-correlated, 
or both vanish. Likewise the Mott lobe at commensu- 
rate bosonic density has vanishing boson superflow and 
nonzero fermion stiffness, is nonzero upon emerg- 
ing from this lobe where the balance between boson- 
boson and boson-fermion repulsions opens a superfluid 
window, with anti-correlated superflow. The superfluid- 
ity between the two Mott regions may be of a novel type 
where the bosons travel along with the fermions (chosen 
to have relatively low density in this work). As a conse- 
quence, the superfluid onset temperature is significantly 
reduced. Finally, we have discussed the signatures of the 
above phases in the momentum distribution function of 
fermions and bosons, which can be measured in time of 
flight experiments. 
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